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ABSTRACT 

Following the recent outburst of the recurrent nova RS Oph on 2006 Feb 12, we measured its near- 
infrared size using the IOTA, Keck, and PTI Interferometers at multiple epochs. The characteristic 
size of ^3 milliarcseconds hardly changed over the first 60 days of the outburst, ruling out currently- 
popular models whereby the near-infrared emission arises from hot gas in the expanding shock. The 
emission was also found to be significantly asymmetric, evidenced by non-zero closure phases detected 
by IOTA. The physical interpretation of these data depend strongly on the adopted distance to 
RS Oph. Our data can be interpreted as the first direct detection of the underlying RS Oph binary, 
lending support to the recent "reborn red giant" models of Hachisu & Kato. However, this result 
hinges on an RS Oph distance of J;540pc, in strong disagreement with the widely-adopted distance of 
~1.6 kpc. At the farther distance, our observations imply instead the existence of a non-expanding, 
dense and ionized circumbinary gaseous disk or reservoir responsible for the bulk of the near-infrared 
emission. Longer-baseline infrared interferometry is uniquely suited to distinguish between these 
models and to ultimately determine the distance, binary orbit, and component masses for RS Oph, 
one of the closest-known (candidate) SNIa progenitor systems. 

Subject headings: techniques:interferometric, stars: novae, infrared:stars, stars: individual (RS Oph) 



1. INTRODUCTION 

Most astronomers are familiar with Classical Novae, 
exploding stars in which an accreting white dwarf (WD) 
in an interacting binary system accumulates enough ma- 
terial for it to become unstable to hydrogen burning. The 
expanding blastwave from one such event Nova Aql 2005 
was recently obse rved by infrared (IR) interferometry 
ijLane et alJl2005f) and the geometric distance was esti- 
mated based on velocities from spectral line observations. 
This result is consistent with the "optically-thick fire- 
ball" model which has been successfully used for twenty 
years to explain the time-evolu tion of the sp ectral energy 
distribution of classical novae lIGehrzlllflSsh . 

While classical novae are expected to recur, very few 
actually have in recorded history. RS Oph is one of 
the handful of so-calle d "recurrent n ovae" with (now) 
6 outbursts since 1898 l)Warnerl fl976). T he most recent 
outbu rst occurred on February 12, 2006 l|Narumi et al.1 
I2006D . and this unusual event motivated intense moni- 
toring by the IR interferometry community. 

The special nature of RS Oph is thought to stem from 
two causes. First, the WD is likely extremely close to 
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the Chandrasekhar limit, since the amount of hydrogen 
needed to trigger an outburst decreases dramatically as 
the WD mass increases. Indeed, detailed models indi- 
cate the WD mass is within 1% of explo ding as a type 
la supernova fe.g.. lHachisu fc Katoll200l(l . Second, the 
mass-losing companion for RS Oph is a red giant (RG) 
with a wind, providing a high density medium for accre- 
tion onto the WD as well as for t h e exp loding blastwave 
to interact with. iBode fc Kahnl l)1985|) have produced 
the most successful model for recurrent novae, drawing a 
clear analogy to extragalactic supernovae and explaining 
t he radio and X-ray light curves in this context. 

lEvans et all l)198S() were the first to study in detail 
the IR time evolution of a recurrent nova. They mon- 
itored closely the 1-3.5 micron flux of RS Oph for about 
3 years after the 1985 eruption. They found the light 
curve had a characteristic (2-mag) decay time scale of 
about 30 days, and comp ared their observation s to the 
generic predictions of the IBode fc Kahnl l)1985fl model. 
They concluded that their observations could come from 
the hot, post-shock gas - this model would predict that 
the IR source should be seen linearly expanding at a 
rate of about ^1 milliarcseconds per day at a distance of 
1.6kpc. This distance estimate is based (most securely) 
on the expanding size of the r adio emission observed i n 
1985 bv lHiellming et alJ (|1986j) and lTavlor et alJ l)1989j) . 
assuming association with the forward sh ock (new ra- 
dio data re-confirm the 1 985 observations; iRupen et alJ 
l2006UO'Brien et alJl2006|L As will become clear, the dis- 
tance to RS Oph is key to the interpretation of the IR 
interferometry data presented here. 

Cha llenging this interpretation, lHachisu fc Katol 
l)2001l see also Kato 1991) have recently produced a 
comprehensive model for recurrent novae meant to 
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explain a wide range of the known novae properties, 
and makes a specific prediction for the origin of the 
near-infra red (NIR) con tinuum which is very different 
from lEvans et all l)1988|) . Following onset of the ther- 
monuclear runaway of the hydrogen shell around the 
WD, the shell expands to AU size, in effect turning 
the WD back into a red giant. The shell stably burns 
hydrogen for a few weeks, shrinking back to the size 
of white dwarf. According to this model, hot post- 
shock gas pla ys no role in forming th e IR continuum. 
Furthermore, lHachisu fe Kato| l)2001|) prefer a much 
closer distance of 600 pc implying a binary separation 
of 2.9 milliarcseconds, easily detectable with current 
interferometers. 

In this Letter, we report first-ever size measurements 
for RS Oph in the NIR using long-baseline interferom- 
etry. Our results are surprising, ruling out the favored 
expanding fireball model, raising doubts about the estab- 
lished distance to RS Oph, and motivating a new model 
for the NIR emission. 

2. OBSERVATIONS 

In the Letter we report on data from 3 different inter- 
ferometers and a summary of observations can be found 
in Tabled Here we briefly introduce each dataset. 

Most of our data were o btained at the Infra red-Optical 
Telescope Array (IOTA, iTraub et all 120031) which has 
baselines between 5-38m. The IONIC3 combiner 
( Bcr ger et all [2003) was used to measure 3 visibilities 
(V 2 ) and 1 closure phase (CP) simultaneously in the 
broadband H band filter (Ao = 1.65^im, AA = 0.3^m). 
Data analysi s procedures have been do cumented in re- 
cent papers ijMonnier et all 12004 l2006j) . For the data 
here, we have adopted a calibration error AV 2 = 0.05 
(relative error) on baselines AB and AC and AV 2 = 0.10 
for baseline BC. 

The Keck Interferometer (KI) was used to observe 
RS Oph on a single night 4 days after the burst (UT 
2006 Feb 16) with a baseline of ~85m. Facility and in- 
strum ent descriptions can be found in recent KI publica- 
tions (Colavita et al. 20 03tlColavita fc Wizinowichl l2003l 
12000^ . The K-band (A = 2.18/Ltm, AA = 0.3/zm) data 
reported here were the by-product of a nulling obser- 
vation (being prepared for a separate publication) and 
the calibration sequence is somewhat modified from the 
standard procedures; in particular, only one ratio mea- 
surement was made per integration. 

Lastly, the Palomar Testbed Interferometer (PTI) ob- 
served RS Oph on 3 nights in 2006 April using the K- 
band detection system (Ao = 2.20/im, AA = 0.4/im) with 
a ~85 m baseline (oriented in NE direction). Detailed 
instrument and data anal ysis descriptions for PTI can 
be fo und in the literature l)Colavitalll999t iColavita et al.1 
1999). Because of the inherent faintness of the source in 
April, coherent integration was used for analysis and a 
large calibration error of AV 2 =0.10 (absolute, not rela- 
tive) was added in quadrature with the internal error for 
model fitting in this paper. 

We have split the data into 3 epochs - from 2006 Feb 
16 to 23 (Days 4-11), from Feb 26 to Mar 13 (Days 14- 
29), and from April 2 to April 18 (Days 49-65). The 
UV-averaged visibility data for each epoch are presented 
in Figure ^ along with some Gaussian profiles for com- 
parison. The IOTA closure phase results are shown in 



Figure [3 also split into the three epochs. 

All V 2 and closure phase data are available from the 
authors; all data products are stored in the FITS-based, 
optical interferometr y data exchange fo rmat (OI-FITS), 
recently described in lPauls et al.1 l|2005lh 

3. ANALYSIS 

Inspection of Figure Q reveals that the visibility curve 
for RS Oph changes very little between day 4 and 
65 (since the outburst). This is surprising since the 
H band brightness faded by a factor >10 during this 
time. This result is discor dant with the generic predic- 
tion of lEvans et alJ l)1988fl who interpreted the IR light 
curve in terms of the time-evolution of post-shock gas at 
10 5 K moving at the speed of the contact discontinuity, 
1400 km/s (following the earlier work bv lBode fc Kahnl 
IT985h . This model requires the IR emission to be seen 
expanding at a rate of ~1.0 milliarcseconds per day (as- 
suming a distance of 1.6kpc), an interpretation that is 
now ruled out. 

Before discussing alternative models in S^l we wish to 
carry-out some model fitting to the interferometry data. 
Here we only consider two simple models - a circularly 
symmetric Gaussian and a binary star model. For all fits 
and calculation of reduced x 2 , we have used the original 
data points before UV-averaging. 

First, we fit a circularly symmetric Gaussian to each 
epoch of data, split by wavelength. Table [3 contains the 
best-fitting Full- Width at Half-Maxima (FWHM) and 
the reduced % 2 for both the V 2 and CP (the model CP 
is always zero for a Gaussian profile). The Gaussian 
model is a reasonable fit for the IOTA visibility alone, 
but clearly cannot fit the non-zero closure phase seen in 
March and April. Also, Figure ^ shows that no good 
fit was possible when combining IOTA with Keck and 
PTI data, indicating this model is too simplistic to ex- 
plain the full range of baselines and/or the wavelength- 
dependence. Sizes derived from the longer-baseline K- 
band data are systematically smaller than those derived 
from shorter-baseline H-band data (IO TA) . 

As discussed in the Introduction, lHachisu fc Katol 
l)2001j) suggest that the nova's IR light curve might be 
due primarily to a rapid increase in brightness of the 
WD as it returns, albeit briefly, to a red giant phase. 
Motivated by this work, we realized that the IR emission 
might be due to the underlying RS Oph binary itself, 
and this might explain the general puzzling features of 
our data: non-expanding emission size, the inadequacies 
of the Gaussian fit, and the non-zero closure phase. 

In order to test this idea, we fit binary models to the 
data for each of the three time periods, treating the 
brightness ratio as independent of wavelength in order 
to fit the H and K band data together. The IOTA, KI, 
and PTI complement each other in Fourier coverage and 
an exhaustive grid search of separations less than 10 mas 
found unique binary star solutions 13 . Table contains 
the best-fit binary models for the three time periods, in- 
cluding the reduced % 2 (V 2 , CP). All three epochs are 
reasonably fit by a similar binary model. The only pa- 
rameter that changed significantly between the epochs 

13 Our limited (u,v) coverage admit some unlikely additional 
solutions with larger binary separations which will be discussed 
fully in a future modeling paper. 
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was the brightness ratio. The closure phase predictions 
for the binary models are plotted along with the clo- 
sure phase data in Figure |2 We note that our uv-plane 
is missing coverage in the NW direction and thus elon- 
gated structure in this direction would be observed fore- 
shortened. 

4. DISCUSSION 

Because the expanding fireball model fails to explain 
the nearly static size scale of the IR emission, we now seek 
suitable alternative emission mechanisms for the time- 
variable IR emission from the recurrent nova RS Oph. 
We have pursued t he reborn red giant (th ermonuclear 
runaway) model of Hachis u fe Katol 1)2001)) and found 
that indeed our 3-interferometer combined dataset can 
be explained by a simple binary model with separation 
of -3.2 milliarcsecond, PA 30° E of N, and a bright- 
ness ratio varying from 2.5:1 to 5:1. Next we subject the 
binary hypothesis to further scrutiny. 

4.1. Binary Interpretation of Near- Interferometry Data 

Bas ed on single-line radial velocity data, iFekel et all 
( 2000) finds the RS Oph binary orbit to be roughly cir- 
cular with a period of 455.72±0.83 days and mass func- 
tion / = 0.221 ± 0. 038 Mo. RGs in symbiotic systems 
are typically 1-3 M© l|Dobrzvcka fe KenvorilT994j) and we 
expect recurrent novae to contain a Chandrasekhar mass 
WD (1.4Mq); these facts combined with the known mass 
function rule out RG masses greater than 2 M© . Assum- 
ing the RS Oph system mass to be 2.4-3.4 M Q , we find 
the component separation to be 1.55-1.74 AU, or (unpro- 
jected) 2.59-2.90 mas a t the 600pc distance preferred by 
IHachisu fc Katd 1)2001)) - only slightly smaller than our 
observed separation of 3.2 mas. 

Since the RS Oph outburst took p lace only 1 mont h 
before maximum redshifted velocity l)Fekel et al.l [20001. 
our measured binary parameters represent the true or- 
bital semi-major axis and orbital f2 for RS Oph with 
only weak sin(i) dependencies. Thus, a small reduction 
in the distance estimate (540 pc) brings the interferom- 
eter binary model in agreement with expectations from 
Kepler's laws. 

The binary model fits (Table |2J) show evidence for a 
change in the brightn ess ratio over time. While the 
IHachisu fc Katd l)2001)) theory predicts a time-changing 
brightness ratio, it is beyond the scope of this paper to 
test the compatibility with the observed IR light curves 
due to complications from the role of the irradiated RG 
photosphere and the presence of a post-outburst WD ac- 
cretion disk. 

4.2. Circumbinary Reservoir of Hot Gas 

The distance estimate of <540 pc derived in the last 
section stands in strong contrast to estimates more 
commonly adopted in the literature. The most sig- 
nificant constraints on distance are set by resolved 
radio observations of the previous and c urrent burst 
lIHieUirrmg et all 119861 iTavlor et alJ 119891 IRupen et all 
2006; IQ'Brien et al]l2006TT By assuming that observed 
radio proper motions (on the sky) can be ascribed to 
the fast-moving ejecta or forward shock, workers consis- 
tently derive a distance of —1600 pc. Similarly large dis- 
tances were found considering interstellar UV absorption 



lines llSniiderslli987D and HI absorption measurements 
(jHiellming et al.111986)) . 

Given the strength of the evidence, we now consider the 
implications of the d = 1600 pc distance. This distance 
would rule out the binary interpretation of the near-IR 
interferometry data laid out in t)4.1l given existing bi- 
nary constraints. Instead, we hypothesize that the IR 
emission arises from a quasi-stationary 14 hot gas reser- 
voir that contributes a combination of emission lines and 
free-free/bound- free emission in the NIR bands. The ob- 
served FWHM of -3 mas is -5 AU at 1600 pc, about 
3x the expected RS Oph binary separation. This size is 
reasonable for a circumbinary disk or reservoir of hot gas, 
perhaps kept ionized by the outward moving blastwave 
or soft X-ray luminosity from the WD itself following 
outburst. This gas reservoir might be analogous to the 
"fallb ack disk" inferred to form after some supernovae 
(e.g.. lWang et~alll200^ . 

Clearly, the hypothesized gas reservoir must be elon- 
gated and somewhat off-center with respect to the cen- 
tral source in order to fit the combined IOTA, KI, and 
PTI interferometry data, especially the non-zero closure 
phases. It is beyond the scope of this Letter to investi- 
gate the details here, and we defer development of this 
model to a future paper. 

4.3. Future Work 

While we have ruled out the important class of ex- 
panding fireball models for explaining the IR emission 
from the recurrent nova RS Oph, more work lies ahead 
to test the other emission mechanisms discussed in this 
Letter. A future study will attempt to synthesize a self- 
consistent model that can explain the time evolution of 
the IR spectrum, NIR and mid-IR interferometer data, 
and multi- wavelength light curves at the same time 15 . 

If the close distance d;$540 pc is confirmed, we have 
a spectacular opportunity to study in detail a likely 
SNIa progenitor and to learn about unexpected shock 
physics controlling the non-thermal radio emission. Al- 
ternatively, the further distance d ~ 1600 pc suggests we 
have discovered a significant and new component to the 
RS Oph Nova remnant and future work will characterize 
the hot circumbinary gas reservoir for the first time. 

We thank J. Sokoloski for productive discussion and 
comments on the manuscript and members of the 
AAVSO IR group for early-time photometry. R.K.B. 
acknowledges NASA and Dr. Bill Danchi for support- 
ing this work. We acknowledge SAO, NSF, NASA, 
CNRS/CNES (France), LAOG, IMEP, and LETI for 
support of IOTA and IONIC3. KI was developed and is 
operated by JPL, MSC and WMKO with funding from 
NASA. PTI was developed by JPL and is operated by the 
MSC on behalf of the PTI collaboration. This research 
made use of SIMBAD, ADS, MSC resources, CHARM2, 
and the 2MASS catalog. Some data presented herein 

14 It is possible to fit our data with a an expanding wind or jet 
component, but this requires fine-tuning the relative proportions 
of multiple components and/or a very asymmetric jetlike emission 
oriented perpendicular to our long (northeast) baselines. These 
possibilties will be investigated more thoroughly in future work 
with an extended dataset. 

15 Note added: future modelin g should address the asymmetric 
radio nebula and jet observed by|U'Bricn et alj 120061) . 
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TABLE 1 
Observing Log for RS Oph 



UcLys since 


Utile 


interferometer 


... , , 

Wavelength 


ZUUO r(!D 1Z 




(Configuration) 


(pm) 


4 


2006Fcbl6 


IOTA (A20B15C00) a 


1.65 


4 


2006Fcbl6 


Keck b 


2.18 


11 


2006Fob23 


IOTA (A20B15C00) 


1.65 


14 


2006Feb26 


IOTA (A35B15C10) 


1.65 


19 


2006Mar03 


IOTA (A35B15C10) 


1.65 


20 


2006Mar04 


IOTA (A35B15C10) 


1.65 


22 


2006Mar06 


IOTA (A35B15C10) 


1.65 


25 


2006Mar09 


IOTA (A35B15C10) 


1.65 


29 


2006Marl3 


IOTA (A35B15C10) 


1.65 


49 


2006Apr02 


Palomar Testbed (NW) C 


2.20 


59 


2006Aprl2 


IOTA (A35B15C10) 


1.65 


60 


2006Aprl3 


IOTA (A35B15C10) 


1.65 


63 


2006Aprl6 


Palomar Testbed (NW) 


2.20 


65 


2006Aprl8 


Palomar Testbed (NW) 


2.20 



Note. — Uniform Disk (UD) diameters of calibrators were generally estimated using getCal, an SED-fitting routine maintained and distributed 
by the Michclson Science Center I http://msc.caltcch.cdul. 



a IOTA used th e following calibrator s: HD 152601 (1.6±0.4 mas), HD 164064 (1.6±0.5 mas), x UMa (3.24±0.04 mas: IBorde et alj|2002l . p Boo 
(3.72±0.12 mas: Ivan Belle et alj[l99gh . HD 143033 (1.9±1.5 mas), HD 156826 (0.6±0.2 mas), HD 157262 (1.5±0.5 mas) 
b KI used the following calibrators: \ UMa (3.35±0.17 mas: ICohen et al!1999l . p Boo (3.92±0.19 mas: ICohen et alj[l99Sh . 
C PTI used the following calibrators: HD 164064 (1.6±0.5 mas), HD 161868 (0.7±0.1 mas). 



were obtained at the Keck Observatory, operated by a scientific partnership among Caltech, UC, and NASA. 
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TABLE 2 
RS Oph Model Fitting Results 






Model 2006 Feb 16 - 23 2006 Feb 26 - Mar 13 
Parameter (Days 4-11) (Days 14 - 29) 


2006 Apr 02 - 18 
(Days 49 - 65) 


Gaussian Profile (fitting only to 1.65//m) 


FWHM (milliarcseconds) 
Reduced x 2 (V 2 ) 
Reduced x 2 (CP) 


3.30 ± 0.09 
0.6 
1.3 


3.47 ± 0.03 
1.3 
3.6 


2.87 ± 0.07 
1.1 
5.4 


Gaussian Profile (fitting only to 2.2/im) 


FWHM (milliarcseconds) 
Reduced x 2 (V 2 ) 
Reduced x 2 (CP) 


2.56 ± 0.24 

ga 

N/A 


N/A 
N/A 
N/A 


2.00 ± 0.09 
1.0 

N/A 


Binary Modcl b (fitting to 1.65/im & 2.2/im data) 


Separation (milliarcseconds) 
Position Angle c (degs E of N) 
Brightness Ratio c 
Reduced x 2 (V 2 ) 
Reduced x 2 (CP) 


3.13±0.12 

36±10 d 
0.42±0.06 
0.6 
1.1 


3.23±0.13 

45±5 
0.40±0.06 

1.7 

1.3 


3.48±0.23 

27±5 
0.21±0.03 
1.2 
0.5 



a Only one Keck data point available for fitting. 

^Brightness ratio assumed the same for H and K bands. Individual components have adopted UD diameters of 0.5 mas. Here, we considered only 
solutions with separations <10 mas. 
c Faintcr component with respect to brighter component. 

^180° ambiguity since closure phase data are indistinguishable from zero for this epoch. 
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Fig. 1. — This figure shows UV-averaged visibility data for RS Oph split into three different time periods. The x-axis shows the spatial 
frequency (projected baseline in units of wavelength) while the y-axis shows the visibility-squared. Four curves representing Gaussian 
profiles arc also included to show the characteristic size and to allow intercomparison of data in the different panels. For reference, all data 
shortward of 25MA derive from IOTA, while longer baseline data come from KI (Epoch I) and PTI (Epoch III). 
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Fig. 2. — This figure shows IOTA closure phase data for RS Oph split into three different time periods. The x-axis shows the hour angle 
of the observations while the y-axis shows the observed closure phase. The solid line shows expected closure phase for the binary model 
parameters found in Table 2 and discussed in J3] The dashed line shows the expected closure phase signal for the symmetric Gaussian 
model. 



